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ABSTRACT 

In cryopumping systems, the ratio of gas condensed to the area 
of the condenser should be made as large as possible. By designing 
an array which would focus the molecules, the condenser area might 
be reduced, and the ratio thereby increased to a satisfactory value. 
The present study was initiated to determine if such focusing could 
be obtained with simple geometric configurations. The kinetics Of 
molecules under free -molecular -flow conditions were examined for 
circles, ellipses, parabolas, hyperbolas, and triangles. This study 
consisted of both an exact theoretical analysis and a Monte Carlo 
computer technique. It was found that the best focusing could be 
obtained with the circular array, which concentrates 50% of the 
molecules into 46% of the area. Since this focusing is so slight, 
an extension of the study is planned to include incident beams which 
enter the array at an angle rather than directly, as in the present 
study. Preliminary results from this extension are encouraging. 
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A STUDY OF CRYOPUMP CONFIGURATIONS IN 
FREE MOLECULAR FLOW REGIONS 

By J. 0. Ballance, W. K. Roberts, and D. W. Xarbell 


SUMMARY 

In cryopumping systems, the ratio of gas condensed to the area of 
the condenser should be made as large as possible. By designing an 
array which would fo^us the molecules, the condenser area might be 
reduced, and the ratio thereby increased to a satisfactory value. The 
present study was initiated to determine if such focusing could be 
obtained with simple geometric configurations. The kinetics of mole- 
cules under free-molecular-flow conditions were examined for circles, 
ellipses, parabolas, hyperbolas, and triangles. This study consisted 
of both an exact theoretical analysis and a Monte Carlo computer tech- 
nique. It was found that the best focusing could be obtained with the 
circular array, which concentrates 50% of the molecules into 46% of the 
area. Since this focusing is so slight, an extension of the study is 
planned to include incident beams which enter the array at an angle 
rather than directly, as in the present study. Preliminary results 
from this extension are encouraging. 


I. INTRODUCTION 

The cryopump has been established as a necessary component of 
vacuum systems in which very low densities and relatively high mass 
flows are required. Since the condensing efficiency of a cryopumping 
system depends more on the physical arrangement of thermal shields than 
on the properties of the gases striking the condensing surface, many 
cryo-array configurations have been developed. The performance of these 
arrays has been examined both theoretically and experimentally. Such 
examinations are usually conducted in the free-molecular-flow region, 
where two basic assumptions can be made: (1) the molecules interact 

only with the surfaces (not with each other), and (2) the molecules 
are reflected from the surface diffusely. Although an exact analysis 
requires knowledge of the accommodation coefficients, sticking coef- 
ficients and capture coefficients, some arrays (such as the chevron 
type) are nearly independent of these parameters. These arrays, however. 
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usually have a very low pumping efficiency; i.e. only a small percentage 
of the molecules entering the array is condensed. 

If there were some way in which the gas molecules could be directed 
or focused, much better pumping characteristics could be obtained for a 
given condensing area. The present study was initiated, using the 
assumptions of free molecular flow, to determine if such focusing could 
be obtained with simple configurations. It was noted that, even if no 
focusing were obtained, the results of this study might be useful in 
determining some of the molecule-surface interaction parameters. 

This report includes an outline of the methods used, a summary of 
the results to date for several simple configurations, a discussion of 
the accuracy of the data, and mention of a possible experimental appli- 
cation of the results. 


II. METHOD 

Typical in any study of free-molecular-flow problems are two 
assumptions: (1) the mean free path of the molecules is much larger 

than any dimensions involved, and (2) molecules colliding with a sur- 
face are reflected according to the Lambert Cosine Law distribution. 

In this study, it was also assumed that the problem could be reduced to 
two dimensions by considering the surfaces to be infinitely long cylin- 
ders. With this assumption, each array was then studied only in cross 
section. A study by a group at Arthur D. Little, Inc. has shown that 
the cosine law can be used in two dimensions without introducing any 
error. 

With these assumptions, two different approaches were used to 
determine exit distributions of molecules reflected from several simple 
geometrical configurations. The first was an exact theoretical approach; 
the second, a Monte Carlo approach, used an IBM 7090 computer. It was 
possible to obtain an exact theoretical expression for the exit distribu- 
tion in only one case, viz., molecules reflected after one collision 
with a circular surface, assuming an incident distribution completely 
uniform over the surface. This calculation is given in Appendix A. For 
other configurations and incident distributions, the final results were 
in the form of integrals which could not be calculated in closed form. 
However, the theoretical approach did yield other important information: 
for the case of the circle, distributions over the surface could be 
calculated as a function of collision number. Expressions for these 
distributions are given in Appendix B and are plotted in Figure 1 for 
two cases: (1) incident flux uniform over the surface (’’random angle”), 

and (2) incident flux uniform across the opening (’’random y”). From 
these distributions, the fraction of molecules exiting the array after 
n collisions, 0^, could be calculated exactly and compared with the 


computer results. This comparison is given in Table I. In this table 
are also given some results obtained by graphically integrating the 
theoretical expression found for the "random y" case. 

TABLE I 

a — Fraction of Particles Exiting Array After n Collisions Circles 


a. Random Angle Incident Beam 


n 

A 

B 

| A-B | 

% 

a 

n 

theory 

a 

n 

computer 

i 

0.6366 

0.6361 

0.0005 

0.1 

2 

0.6259 

0.6295 

0.0036 

0.6 

3 

0.6279 

0.6187 

0.0092 

1.5 

4 

0.6281 

0.6012 

0.0269 

4.3 


b. Random y Incident Beam 


n 

A 

B 

C 1 

H 

% 

H 

% 

a 

n 

theory 

a 

n 

computer 

a 

n 

graphical 

a 

0.6667 


0.6756 

0.0059 

0.9 

0.0089 

m 

— 

0.6186 


0.6757 

0.0019 

0.3 

0.0571 

9.2 

Hi 

0.6298 

0.6400 

0.6759 

0.0102 

mm 

0.0461 

m 

D 

0.6285 

0.6538 

0.6774 

0.0253 

4.0 

0.0489 

7.8 


It is seen that the computer results for the fraction are in error 
by a maximum of 4.3% while those obtained by the graphical method are 
in error by as much as 9.2%. The graphical method was therefore dis- 
carded in favor of the computer method. 

The computer method, using the Monte Carlo technique, consisted of 
the following steps: 

(1) A random number, R ly between 0 and 1, was generated by the 
computer . 
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(2) From R 1? a random number y L = l-2R x (between -1 and 4-1) was 
computed and used for the y coordinate of the incident particle. 

jj3) Using the general equation for the geometry of interest, the 
x coordinate, Xi, was then calculated. 

(4) The equation of the normal to the surface at (x x , y x ) was then 
computed. 

(5) A new fandom number, (also between 0 and 1) was generated 

and an angle, calculated from 6 = sin"-*- (I- 2 R 2 ). 

(6) The equation for the particle line, at an angle 5 with the 
normal, was then computed; this yields the cosine distribution for 
the reflected particles. 

(7) The intersection of this particle line with the y-axis, "Y", 
was then calculated. This value, Y, was then tested and placed in the 
correct group. 

(8) If the particle did not fall into any of the 20 groups from 
Y = -1.0 to Y » 4-1,0, the intersection of the particle line with the 
array was calculated. This yielded a new (x x , y x ) for the particle’s 
second collision and the computer then returned to step (4). 

(9) The program was continued until a particle had either been 
placed in one of the exit groups or had made four collisions (six 
collisions in the case of the triangle), whereupon it was recorded as 
"out" and the program started over at step (1). 


III. DISCUSSION OF RESULTS 

Computer results were printed out in twenty groups, from y = -1.0 
to y - 4*1.0 in intervals of 0.1, with each group broken down into first, 
second, third and fourth collisions. These distributions, for the five 
geometries considered, are presented in Figures 2 through 6, with only 
Figure 2 including the distribution as a function of collision number. 

It is clear that very little focusing is obtained with any of these 
arrays. A measure of focusing is given by the percentage area into 
which the middle 507, of the particles fall. These percentages are 
given in Table II, for each of the arrays considered. 



5 


TABLE II 

Percentage Area into which Middle 50% of Particles Fall 
for Various Surface Geometries 


Geometry 

% Area 

Circle, Random Angle 

51.3 

tt u y 

45.8 

Ellipse, A = 0.5 

47.5 

" , A = 1.25 

46.9 

” , A = 1.50 

46.7 

" , A = 1.75 

46.2 

" , A = 2.00 

46.5 

Parabola, P = 0.5 

47.7 

■ 

H 

46.8 

■H 

48.4 


Geometry 

% Area 

Parabola, P = 0.0833 

47.7 

Hyperbola, A = -1 

48.3 

" , A = -2 

46.9 

" , A = -3 

46.6 

" , A « -5 

4745 

" , A = -7 

47.8 

Triangle, 10° 

48.9 

" , 30° 

48.6 

" , 45° 

47.6 

" , 60° 

47.7 


It can be seen that the best focusing is obtained with the circle 
with random y incident beam, which "concentrates 11 50% of the particles 
into 45.8% of the area. 

Since these results were rather uninteresting, it was then decided 
to examine cases where the incident beam enters the array at an angle p 
with the opening. If better results could be obtained under these condi- 
tions, then several simple arrays might be combined (in a clover-leaf 
pattern, for example) to focus the molecules. Preliminary results of 
this study are presented in Figure 7, which shows the distribution of 
particles after one collision with a circular array, where the incident 
beam is at an angle p with the y-axis. These results are somewhat more 
encouraging, and these studies are being continued at the present time. 


IV. DISCUSSION OF ERRORS 

Two methods have been used to determine the random errors inherent 
in the computer approach. The first of these is based on the asymmetry 
of the results, as follows: Since both the incident beam and the arrays 
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considered were perfectly symmetrical with respect to the x-axis, it was 
expected that the exit distributions across the y-axis would also exhibit 
this symmetry. For example, the same number of particles should fall in- 
to group 20 as into group 1, the same into group 11 as into group 10, etc. 
With the finite number of points (10,000) considered by the computer, there 
was, of course, some asymmetry. Theoretically, the difference between two 
values obtained should be 0 ± 2a (a = standard deviation), with 95% confi- 
dence, or 0 ± .6745a with 50% confidence, where a standard deviation in 
this case is given by ^^Npq’ (N is the number of points considered = 10,000; 
p is the probability of a point falling into a group * ^ qqq ; n is the 
average number in a group; and q = 1-p). Both the percentage probable 
error and two standard deviations are plotted in Figure 8a, where the 
solid curves have been plotted from the above theoretical expressions, and 
the points have been calculated from the computer data. The figure indi- 
cates that the asymmetry obtained was approximately what was expected 
theoretically. 

The second method used to determine the errors in the computer method 
was based on a comparison with the only case solved exactly by theoretical 
analysis: the case of the distribution of particles exiting the array 

after one collision with a circle, where the incident particles are dis- 
tributed uniformly over the surface. Twenty-four different computer runs 
were made, varying the number of incident particles. The probable error 

was then determined from the equation: P.E. = .6745 (Sx?) / 20 ,* where x^ 

is the difference between the number of particles in the ith group and 
the average; the index i runs from 1 to 20, since there are 20 groups. 

The prob able error was also calculated theoretically, using P.E. = .6745a 
and a = a / Npq 1 where N is the number of points considered, p = .03183 (calcu- 
lated from theory) and q = 1-p. The results are shown in Figure 8b. It is 
seen that the computer errors are within those predicted theoretically, and 
that the probable error for 10,000 points is approximately 3.7%. Thus the 
distributions shown in Figures 2 through 6 can be considered accurate within 
2 to 5%, depending upon the number of particles in a group (for 250 particles 
the error is 4.2% or 11 particles; for 600 particles the error is 2.7 % or 
16 particles) . 


V. EXPERIMENTAL APPLICATION: MEASUREMENT OF ACCOMMODATION COEFFICIENT 

As stated earlier, an exact analysis of an array requires knowledge 
of the accommodation coefficients, sticking coefficients, etc. At present, 
this study has not progressed far enough to be concerned with these para- 
meters; however, the results might be very useful in their determination. 

In particular, an immediate application of the results would be the 
determination of the accommodation coefficient. 
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The accommodation coefficient, a, is defined by: 


a = 


T - T 
o n 


T 


o 



where T , T , and T are the absolute temperatures of the incident gas, 
the ref?ected gas, and the reflecting surface, respectively. 

Thus the temperature of Ni molecules after one collision is 

T = T -a(T -T ). It is easily shown that the temperature of N 
l 0 o w n 

molecules after the nth collision is: 

T = T + (T - T ) (1 - a) n (1) 

n w o w 


For a given experimental array, some molecules will exit the array 
after one collision, others will collide more than once with the sur- 
face of interest, and exit after n collisions. Considering all exiting 
molecules, the actual measured temperature of the exiting gas will be a 
weighted average of all these, given by: 



( 2 ) 


where N is the number of molecules exiting the array after n collisions 
n 

The summations are shown here to extend to infinity, but if a cut- 
off is imposed, this can be replaced by a finite number. For example, 
if the cutoff is set at 99.9% (i.e. , 99.9% of the particles have exited) 
the upper limit on the summation is a number between 2 and 50, depending 
upon the geometry. 

The fraction of particles exiting the array after n collisions, Q^, 
has been discussed above (See Section II). 

For n » 1: <*x = N x /No or N x = a x N 0 . 


Similarly, 


N 2 = a 2 (No-Ni) , etc. 


N 

n 




In general. 


(3) 
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Substituting (3) and (1) into (2): 



( 4 ) 


Now using the results of Section II above for a n , assuming a con- 
stant incident gas temperature, T Q , and a constant surface temperature, 
T w , the exit temperature, T can be plotted as a function of the 

accommodation coefficient, a. Then, by experimentally measuring T exit , 
the accommodation coefficient can be found graphically* 

Unfortunately, the present study has shown that the distribution of 
exiting particles across the opening of the array is not uniform and, 
more importantly, the distribution changes with collision number. Thus, 
the measured value of T exit will change as the probe is moved along the 
y-axis (the opening). However, by using the computer data, the various 
°ti can be determined for any given section of the opening and the above 
procedure used to yield the accommodation coefficient. 

Plans are now in progress to conduct this experiment in the near 
future . 


VI. CONCLUSIONS 

This study has proven the great value of the Monte Carlo computer 
technique in cases where exact theoretical techniques cannot be found. 
This technique is accurate to better than 47o when 10,000 points are 
used. 


The study has shown that very little focusing can be expected from 
simple geometric surfaces if the incident beam of molecules is parallel 
to the axis of symmetry. Some focusing is possible when the incident 
beam makes a large angle with the axis of symmetry. 

An elementary application of the results of the study might be a 
determination of the accommodation coefficient by a simple experiment. 
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APPENDIX A 


Distribution across Y-Axis of Particles Emitted after One Collision 
Incident Distribution Uniform over the Reflecting Surface 

("Random Angle 11 ) 


List of Symbols (See Figure A) 


y 


distance along the opening, measured from the center 
of the semicircle. 


P 

« 


point on semicircle at angle Oi measured from negative 
y-axis. 


a 


N 


o 


angle measured from negative y-axis to normal at 
point P . 

a 

total number of incident particles. 


9 

N 

a,cp 


number of particles per radian reflected at point P . 

a 

angle between path of reflected particle and normal 

to surface at P . 

a 

number of particles per radian reflected at point P 
which are reflected at an angle cp. ^ 


0 


angle from normal at P to line from P to the point 

(0,y). a 


N 

y, a 


number of particles per radian from P^ that pass 
through a section from (0,0) to (0,y). 


N 

y 

7 
b 

Consider an incident distribution of N Q particles which is uniform 
over a semicircular surface as shown in Figure A. Assuming that these 
particles are diffusely reflected according to Lambert's Cosine Law, 
the distribution of reflected particles across the y-axis from y = -1 
to y = +1 can be calculated as follows: 


total number of particles through the section from 

( 0 , 0 ) to ( 0 ,y) from all P . 

a 

angle from positive y-axis to normal at P (7 - it - a) . 
distance from P to point ( 0 ,y). 
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y 



Let = total number of particles 
per radian reflected at point P . 


Then N = 


N da 
a 


= itN 


a 


The number of particles per radian reflected at an angle cp from the 
normal to the surface at is given by: 


N = k N cos cp (Lambert's Cosine Law) 
Q ly Cp 01 


To find the constant k: 


N = 

a , 


+ 1 
p / 


jr 

2 


+ f 

2 


hence 


k = — , so that 


N = — N cos cp. 

a, cp 2 a 


Now the total number of particles per radian from p a that passes 
through a section of the opening from 0 to y is given by: 


e e 

*y,a m f 008 * ^ ’ 


2 “a 8ln 6 
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and the total number of particles passing through this section from all 

? a U: 


n 71 

N = / N da = / i N 

y J y, a J 2 a 

0 0 


sin 0 da 


(i) 


Thu’s to calculate N , sin 0 must be written in terms of a. Referring 
to Figure A: ^ 


sin 0 _ sin 7 
y " b 


. ^ y s in 7 

or sm 0 = * — 

b 


but 


and 


sin 7 = sin (ic-a) = sin a 

% 


■0 


b = I 1 + - 2y cos 7 


') 


Also 


cos 7 = cos (it-a) = -cos a 


so that 


b = I 1 + + 2y cos a 


* 


Finally 


sin 0 = 


y sin a 


0 


. \* 

1 4* y^ + 2y cos CM 


Inserting this expression into equation (1): 


i in a da 


n = r 4 n 2-2i> 

y J 2 a _ v 

0 ( 1 + y 2 + 2y cos aj 


This integrates immediately to; 


N / v % -i 

N y = J (^ + / + 2y cos aj 


h ~| n 
0 
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N 

y 


!« r 

2 


(l + y 1 


- 2y) % - 


(1 + y 2 + 2y) 


N 

y 


N r 

a 


±(1 - y) - ± (1 + y) 


Now (1 + y + 2y cos 


a) % = b, 


which is a positive number (a distance). 


Hence , 


N 

y 



(i - y) - (i + y) 



Thus the 
0 to y is 
y-axis is 


total number of particles that pass through any section from 
simply proportional to y; i.e., the distribution across the 
uniform. 
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APPENDIX B 

Distributions of Particles over a Semicircular Surface 
as a Function of Collision Number 


In general, the number of particles per radian making the nth 
collision with a semi-circular surface, at an angle a measured from 
the y-axis, is given by: 


N (n) = - £- N (n) 

a dp w p 


J p=a 


where N^ n) = 2 F (n) (p,p) - F (n) (0,p) - F (n) (rt,p) + N^ n_1) 
P P 


and (a,p) = N^ n ^ (^cos j cos & + sin j sin da 


can be calculated from the previous expression for N^ n 


Case I. Incident particles uniform over the semicircular surface 
(random angle incidence), N = total number of incident 
particles, ° 



N h x 
o 
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H a 5) - (N 0 /128«) [l28 - (93 + ^ - 29a + 4*a - 2 §S . JS£ + 

cos 2 + (-93 + 2Z£ . |f + . 29a . 3 *a + M . *2! + |!) sln 2 

Case II. Incident particles uniform over the opening of the array 
(random y incidence). 



o'lQ 
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Figure i. Distribution of particles over a semi-circular surface 
at first five collisions (cosine law of reflection) 




Figure 2. Distribution of particles across opening of CIRCLES as a function of 
collision number (10,000 incident particles) 














Figure 4. DISTRIBUTION OF PARTICLES ACROSS OPENINGS OF PARABOLAS 























Figure 6. DISTRIBUTION OF PARTICLES ACROSS OPENINGS OF TRIANGLES 
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